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ABSTRACT 


The  mechanism  of  the  doublet  at  supersonic  speeds  is  described  using 
elementary  physical  reasoning  insofar  as  is  possible.  The  cause  of  the 
infinities  introduced  into  the  equations  by  the  differentiation  process 
across  the  Mach  cone  is  discussed.  In  addition  the  physical  process  involved 
in  passing  from  the  fixed  to  the  moving  doublet  is  explained.  Finally  some 
applications  of  doublets  to  supersonic  aerodynamic  problems  are  included. 
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INTRODUCTION 


The  use  of  sources  and  doublets  for  the  analysis  of  flow  fields  has 
led  to  many  useful  results  in  aerodynamics.  In  the  developments  concerning 
incompressible  flow,  the  physical  pictures  of  sources  and  doublets  were 
retained  while  for  supersonic  flow  the  approach  has  been  more  mathematical 
in  the  sense  that  the  source  and  doublet  flows  have  been  referred  to  as 
particular  solutions  of  the  differential  equation  involved  rather  than 
physical  flows.  This  mathematical  approach  has  the  disadvantage  that  the 
development  of  intdi-tion  by  designers  and  engineers  regarding  supersonic 
flows  is  impeded. 

In  reference  1,  the  development  of  supersonic  theory  utilizing  the 
source  i.e.,  a point  in  the  fluid  at  which  fluid  is  added  or  withdrawn  at 
a given  velocity.  The  effect  of  moving  the  source  is  then  derived  following 
the  superposition  method  of  Prandtl  given  in  reference  2.  This  process 
enables  one  to  follow  the  physical  mechanism  involved  in  the  fluid  motion 
around  bodies  moving  at  speeds  faster  than  sound.  Another  concept  des- 
cribed in  reference  I is  the  pressure  source  in  which  fluid  is  introduced 
with  a certain  acceleration.  For  many  problems,  the  pressure  source  seems 
to  give  a simpler  picture  of  the  flow  than  the  velocity  source. 

Some  problems  in  fluid  mechanics  can  be  readily  solved  with  the  use  of 
the  double  source  or  doublet.  In  this  case  a source  and  sinlc  are  brought 
close  to  each  other  but  the  strength  times  the  distance  is  held  constant. 
Mathematically  this  situation  can  be  represented  as  the  derivative  of  the 
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source  along  the  line  joining  the  source  and  sink.  The  pressure  doublet  then 
corresponds  to  the  derivative  of  the  pressure  source.  The  stream  lines  can 
be  readily  visualized  i.  e.:  one  half  of  the  fluid  goes  directly  from  the  source 
to  the  sink  and  the  remainder  flows  from  the  opposite  side  of  the  source  to 
the  corresponding  side  of  the  sink. 

It  has  been  found  in  supersonic  flows  that  the  pressure  of  both  the  source 
and  doublet  is  infinite  along  the  Mach  cone.  This  situation  causes  little 
trouble  in  the  case  of  the  source  but  has  required  special  mathematical  tech- 
nique in  the  case  of  the  doublet.  Another  characteristic  of  the  doublet  is 
the  difficulty  of  passing  from  the  fixed  to  the  moving  doublet..  It  is  shown 
that  Prandtl’s  approach  for  the  source  as  modified  in  reference  1 does  not 
seem  to  be  applicable  to  the  doublet.  One  of  the  purposes  of  this  report  is 
to  clarify  the  physics  involved  in  the  two  difficulties  listed  above.  Another 
purpose  is  to  present  some  of  the  known  applications  of  the  doublet  without 
the  use  of  the  special  mathematical  techniaues  normally  applied. 
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PART  I 


SINGULARITY  AT  THE  MACH  CONE; 

The  singularity  at  the  Mach  cone  for  the  doublet  is  clearly  tied  up 
with  the  fact  that  one  takes  derivative  of  the  source  pressure  which  jumps 
from  zero  to  infinity  on  the  Mach  cone.  From  reference  1,  we  have  for  the 
pressure  of  the  moving  source  the  following  expression: 


e fcr) 


1,1  f ' 4TT|>-(M*iXy%-rz9]  ^ 


where  p = (pressure)  lbs/sq.ft. 

£ = strength  of  source  ( lbs. /ft.) 
x,  y»  z *■  coordinates  in  the  moving  system  (ft.) 

M - Mach  number  of  source  moving  along  thte  x axis  in  the 
negative  direction 

* starting  time  of  source  striking  the  point  at,  y,  z,  at 
time  t 

T=  t - +zx 

c 

For  the  case  where  the  source  strength  is  steady  i.e.,  f(r)  = 1 , 
there  will  be  two  waves  striking  the  point  x,  y,  z at  each  value  of  time 
i.e.,  the  backward  moving  part  of  one  wave  and  the  forward  moving  part  of 
another  wave  emitted  at  an  earlier  time.'  The  steady  state  expression  is 
therefore: 


1.2 
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This  wave  will  accelerate  the  air  particles  according  to  Newton's  Law 
i.e#,  proportional  to  the  pressure  gradient.  Application  of  this  formula 
demonstrates  the  cause  of  the  mathematical  difficulties  in  handling  doublets 
at  supersonic  speeds  i.e.: 


ZL 


Considering  this  expression,  we  see  that 


always  has  the  same  sign 


behind  the  Mach  cone  and  all  particles  would  be  accelerated  down.  However 
as  the  particle  enters  the  Mach  cone  it  vail  be  accelerated  up  since  the 
pressure  jumps  from  zero  to  infinity.  This  upward  acceleration  which  is 
neglected  in  the  usual  computation  cancels  the  infinite  downward  acceleration 
which  gives  rise  to  the  infinite  velocity  obtained  in  the  mathematical  pro- 


cedure as  follows: 


w = - f- 

ATTpu  J 


[xa-CM*-o(yvz»; 


where  the  * refers  to  the  value  of  x which  will  cancel  the  denominator  i.e., 
the  point  where  the  particle  is  at  the  Mach  cone.  The  integral  yields  the 
expression: 


w = 


I ~CS\Z?)  7 


so  that  substitution  of  the  limit  x*  gives  an  infinite  velocity  everywhere. 
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From,  simple  acoustic  theory,  the  velocity  due  to  a plane  wave  is  t/pc. 
Since  the  Mach  cone  appears  plane  to  a particle  very  close,  this  formula 
can  be  used  to  evaluate  the  axially  symmetric  velocity  caused  by  the  positive 
pressure  jump  just  ahead  of  the  Mach  cone.  The  vertical  velocity  caused  by 
this  jump  is: 


1.6  w*  t ZC056  — SsSfl-  L 

p«-  Vv)H2“-  zvp  (o  v'xmms.x^+z*) 

where  ^ is  the  Mach  angle. 

Now  and  =£==*-  ~ so  that  the  infinite 

c oT~  V/ 

velocity  due  to  the  acceleration  just  behind  the  Mach  cone  is  cancelled  out 
by  the  infinite  velocity  due  to  the  pressure  jump  just  in  front. 

The  above  derivation  gives  a physical  explanation  of  the  infinities 
which  arise  in  supersonic  theory  and  which  have  reouired  various  matae.  at i cal 
schemes  such  as  Kaiamard's  method. 

The  velocity  caused  by  a pressure  source  at  the  origin  is  therefore 
given  by  the  expression: 


1.7 


vy 


> e 


xz 


«2.Trpu  iy*+zz)  }/x~  Cm  ~iXy  N- z*) 


The  velocity  caused  by  a doublet  is  .obtained  by  taking  the  derivative  i.e., 

V/  = -£! 

1.8 
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We  can  obtain  emotion  1.8  directly  by  considering  tbi  c.  .:v„  ..aj,  for  the 


pressure  doublet. 


to  _ € 3/ 

r 2TT  [1^-- 


The  vertical  ::ccelerution  is: 


3L 


- ^ o U — - ■[  ^ 1 3/  _l  3(M*.Q*Z2  ) 

~ T dx  ~ iir  ((x^ (ki^X^\z:x)3  a^(xa-(^*-iXyl+z?)  j 


Integration  vdth  respect  to  x yields  for  the  velocity: 


W “ 


2TTpU 


1 LsW; 


5Z*  x _ _L  X3 s ] 

(v.)VZ*)x  x»  3 (xl>Civ\^Vz^ 


If  we  neglect  the  lover  limit,  the  result  is: 


w~ 


2.TTPU  l \/ 


M+Z>)x  tfx*  (M1-  iX'dt2lXlSh-£1X*  - (MViX^Vz4))3^ 

Comparison  of  1.12  with  1.8  shows  that  the  first  term  of  1.12  eouals  the 
second  term  of  1.8.  Combining  the  other  terms  of  1.3  yields: 

1V  3 fa*-*)  Z^HlI 

( X1-  - ( ^x-  ilhN-zh)3^  bf+zr) 1 


3 71 


WADC  1R-52-290 


4 


An  identical  result  is  obtained  from  1.12  which  shows  the  identity  of 
the  two  methods  but  does  not  justify  omitting  the  terms  with  X*  as  the  lower 
limit.  This  justification  is  obvious  however  from  equation  1.6,  the  deriv- 
ative of  which  with  respect  to  21  exactly  cancelling  out  the  infinite  terms 
caused  by  the  lower  limit  in  I. 11. 
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PART  II 


THE  MOVING  DOUBLET 


In  reference  1,  the  effects  of  motion  on  the  pressure  or  velocity  source 
is  demonstrated  in  a very  simple  manner.  A moving  source  is  represented  as  a 
series  of  fixed  sources  emitting  fluid  for  short  intervals  of  time  in  successive 
positions  corresponding  to  the  location  of  the  moving  source.  The  pressure 
change  here  is  caused  by  overlapping  of  the  waves  at  a given  point  in  space. 

This  treatment  essentially  follows  that  of  PrandtT  given  in  reference  2. 

The  usual  treatment  which  yields  the  moving  doublet  is  to  take  the  de- 
rivative of  the  moving  source  although  it  is  possible  that  a clearer  picture 
could  be  obtained  by  directly  passing  from  the  fixed  doublet  to  the  moving 
doublet.  An  attempt  to  use  the  Prandtl  method  for  this  purpose  proved  un- 
successful as  explained  in  the  derivations  of  this  part  of  the  report. 

We  consider  first  the  potential  caused  by  a fixed  source. 


II. 1 


(K  - € f(  t ~rO 

r 4irr 


The  fixed  doublet  then  follows  from  the  expression  for  the  derivative 


II. 2 


gi  r ((t- r/c~) z , f* (tzw. 
4irt  r3  ' ric 


This  is  the  same  expression  as  the  velocity  of  fluid  from  the  source.  The 
first  term  represents  the  velocity  one  v;ould  obtain  from  an  incompressible 
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source.  The  second  term  represents  the  velocity  due  to  the  acoustic  wave 
which  is  emitted  whenever  the  source  changes  its  strength.  For  steady  flow 
behind  the  wave  front,  we  can  neglect  the  second  term.  It  was  noted  in  Part 
I that  the  velocity  <s»nsed  by  this  second  term  cancels  the  infinity  occuring 
in  the  mathematics  for  velocity  calculation. 

If  we  calculate  the  potential  of  a moving  source  using  Prandtl's  method, 
we  obtain  the  following  relation  (see  reference  1): 


H.3 


6 

^rrr 


C r 


fCt-Vc) 


Where  the  term  on  the  right  hand'side  between  the  bars  represents  the  contraction 
of  a wavelet  due  to  its  motion.  In  the  fixed  system  the  wave  is  contained 
between  two  spheres  having  the  same  center  while  in  the  moving  system  the 
center  of  the  two  spheres  is  different.  The  potential  increases  in  proportion 
to  the  contraction  ratio.  We  can  visualize  this  as  follows:  each  wavelets  is 

made  up  of  a number  of  elementary  wavelets  each  of  which  represents  a certain 
value  of  the  potential.  When  the  wavelet  is  contracted  as  a result  of  motion, 
the  elementary  wavelets  overlap  so  the  potentials  add. 

The  value  in  11,3  represents  the  position  of  the  elementary  wavelets 
striking  the  point  x,  y,  z,  at  time  t.  Now  if  we  express  "£|in  terms  of  x,  y, 
z,  and  t.iA.,  - U(X-Vti)  = C fa-Vt) *•- yW)  the 

expression  for  the  potential  in  the  moving  system  X-  x-vt  is: 


II. 4 


€ 

4TT  i)fc»VZ*) 
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which  is  the  veil  known  result. 

When  ve  take  the  derivative  of  with  respect  to  21  , the  term  outside 
the  brackets  in  II. 3 turns  out  to  be  the  velocity  of  the  fixed  source  however 
inside  the  brackets  both  T*  and  t|are  functions  of  Z+.  so  that  these  terms 
must  be  included  in  the  derivative  expression.  This  means  that  the  velocities 
add  because  of  the  overlap  as  is  the  case  for  the  potential  but  additional 
gradients  in  the  potential  exists  bee  use  of  the  contraction  and  because  of 
the  change  in  £|With  £ v.iiich  result  in  additional  velocities.  This  means 
that  we  c-nnot  add  the  velocities  at  a point  caused  by  a succession  of  fixed 
sources  and  obtain  the  velocity  of  the  moving  source  but  ve  must  also  add  the 
incremental  velocities  caused  by  the  motion.  A similar  condition  exists  in 
the  case  of  the  pressures  due  to  a moving  doublet. 

If  we  su  stitute  C(t'ti^for  V-  which  is  its  ecuivalent  we  obtain  for 


from  II. 3: 


II. 5 


C£ 


T - uO-uii) 

The  derivative  then  yields: 

II.6  - - C 6 (u*-C*)  ft 

^ G-kt-tO  -O(x-uto)a 

Using  the  relation  mentioned  above  i.e.: 


H.7  (x~ 


which  gives  after  differentiation: 


WA.'V’  TK-5P-P90 


3 


ii  .a 


at  i 

dz 


so  that: 


II. 9 


d Z 


Z, 

- u(X-uti) 

£ Z (m-0 

dir  (y-vt)*-  <M»,X91+ZJ)) 3/4 


whert  in  the  moving  system  of  coordinates. 

This  part  of  the  report  indicates  that  the  superposition  method  of 
Prandtl's  does  not  seem  to  be  applicable  for  the  velocity  due  to  a source  or 
the  pressure  due  to  a doublet.  However  there  is  no  question  that  a moving 
doublet  is  the  equivalent  of  a series  of  fixed  douolets  each  emitting  suc- 
cessive pulses  from  the  position  of  the  moving  doublet.  The  problem  here 
seems  to  obtain  the  correct  value  for  the  summation  of  -the  velocities  and 
pressure  pulses.  For  a more  detailed  understanding  of  the  details  of  such 
a summation,  it  appears  that  more  research  is  needed. 
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PART  III 


APPLICATIONS  OF  DOUBLETS  TO  SUPERSONIC  PROBLEMS 


The  first  step  in  the  application  of  the  doublet  to  the  solution  of 
aerodynamic  problems  at  supersonic  speeds  is  to  relate  the  doublet  strength 
to  a quantity  at  the  surface  of  the  airfoil.  It  will  be  shown  by  analogy 
with  the  velocity  source  that  the  strength  of  a pressure  doublet  is  eoual  to 
the  lift  on  an  element  of  the  wing  surface.  The  pressure  of  a doublet 

corresponds  to  the  velocity  of  source  since  the  first  is  obtained  by  the 
partial  derivative  of  the  pressure  source  and  the  second  by  the  partial 
derivative  of  an  ordinary  fluid  source.  In  a fluid  source  the  potential  and 
velocity  are  riven  respectively  by: 


III.l 


III. 2 


& = — - — 

_ _ c z_ 

az.  ~ r 


Equation  III.l  and  III. 2 correspond:  to  a sink  since  the  vertical 

velocity  is  negative  for  positive  X-  anh  positive  for  negative  X • order 
to  evaluate  £ , x^e  assume  a small  area  is  represented  by  a distribut- 

ion of  sinks  of  uniform  strength  per  unit  area  i.e.,  £ is  replaced  by 
Very  close  to  the  surface  the  flo\v  will  be  plane  so  that  the  velocity  is  given 
by  the  volume  flow  per  second  divided  b,  the  area  i.e. 


III. 3 


w 


cfx  A* dy 
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Since  W corresponds  to  the  pressure  £ is  the  twice  the  pressure  on 
the  upper  surface  and  £ is  the  lift  acting  on  the  element  ^3.  In 

corresponding  fashion  6 for  a velocity  doublet  can  be  shown  to  be  equal 
to  the  difference  in  potential  times  the  area. 

From  the  above,  we  have  a relation  between  the  lift  on  an  element 
on  the  surface  of  a wing  and  the  pressure  in  the  field.  Using  Newtons  Law 
we  can  compute  the  acceleration  of  the  fluid  from  the  pressure  gradient  and 
compute  the  velocity  by  integration  with  respect  to  time. 

We  will  consider  as  the  first  example  of  the  method,  the  derivation  of 
the  formulas  for  the  downwash  due  to  a lifting  doublet  at  both  subsonic  and 
supersonic  speeds.  The  pressure  caused  by  a doublet  is  given  by  the  expression 
derived  in  Part  I with  substituted  for  £ i.e. 

+ UMQ  dj Jn.  F Xr)  d7"/d  z 

Where  L<H<1  h=lift  on  the  element  JjJn  (lbs) 

f(r)  =•  function  showing  lift  variation  with  time 
"T**  •=  time  of  emission  of  pressure  pulse  striking  the  point 
x,  y,  z at  time  t. 

Considering  first  the  steady  state  case  f(r)  is  the  constant  unity  and  7”*  is 
given  by  the  relation  (Reference  X.) 

in.  5 c\t - rf  = w (t-kSf-f  ( z3" 
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The  left  hand  side  of  IIL5  equals  the  distance  the  wave  travels  in  the 
elapsed  time  (t-T)  and  trie  right  side  shows  the  space  distance  from  S,n. 
to  the  point  x,  y,  z.  Solution  of  this  equation  shows  only  one  real  root  for 
subsonic  speeds  but  two  for  supersonic  speeds.  This  means  that  two  pressure  out 
waves  reach  each  point  at  supersonic  speeds  where  as  only  one  occurs  at  sub- 
sonic speeds.  So  m*  * for  m>1  and  for 

Newtons  Law  gives  O X—  ~ or  for  the  steady  case  x ij 

r dt  ai  At  6% 

Writing  this  equation  for  the  subsonic  case  give  for  W > the  following  expres- 
sion at  1—0  the  plane  of  the  airfoil. 


III. 6 W (X)'j)  - 


/■% 

La$<U  (Mh)  f djc 

4-itpu  1 (of-S)'-  cm -iXv-nr-^ 


and  after  integration 


III. 7 V/  = 


L <JS<W  f L 


Far  behind  the  doublet  we  obtain: 


III. 8 


Ljjjn  r > 

airpu  L es-*)1 


An  if  we  integrate  this  equation  for  — Y*  to +**  we  note  immediately  the 
similarity  between  a row  of  doublets  and  a horsehoe  vortex  extending  from 

“ YV  to  4 ft  . 

For  the  supersonic  case  we  use  2 for  the  value  and  use  the  Mach 
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line  as  the  lower  limit.  This  yields: 


III.9  W(x,vj) 


L djtln  f (m -Q  dx 

iTfpu  Jx*  -( MViXv-n)^3- 


It  is  shovm  in  Part  I,  that  the  lower  limit  of  this  integral  gives  a 
downwash  of  zero  so  that: 


III. 10 


W(*,^)  = 


U dii  in  1 

| CK-0  1 

a.Trpu 

L^J^s  y-  - 

If  ) is  large  compared  to  ( S"*),  III. 10  ecuals  III .8,  showing 

the  downwash  at  far  distances  from  both  a supersonic  and  subsonic  doublet  is 
the  same. 

The  next  problem  which  will  be  considered  here  is  the  conically  loaded 
delta  wing.  It  is  convenient  to  use  for  this  problem  the  vertical  acceleration 
due  to  a line  of  doublets  rather  than  the  vertical  velocity.  We  obtain  from 
formula  III. 9,  the  result 


p w _ L ( M ~0 3/ 

dx  aTTTpU  1 


We  use  the  coordinates:  “K  ecual  the  tangent  of  the  angle  from  the  center- 
line  and  ^ eeual  the  distance  to  the  point  along  the  ^ axis.  The  vertex  of 
the  delta  wing  is  at  the  origin  and  its  centerline  is  on  the  )C  axis.  Therefore 

so  that: 

»* 

. LOflah,  P fM'-Olcii 

Jo  (&-IT- 2 


n ~ 

and 

III. 12 

d V/ 
d* 
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The  integral  here  has  the  form 


$di/x^x 


X - a + cf*- 


This  yields: 


III. 13 


_ _ UhHfc  f 

2>Y  ~ ZTTp  U / ) 


The  infinity  along  the  Mach  Cone  of  the  last  doublet  to  influence  a*: 

&TC 

is  thrown  out  since  it  was  proved  in  Fart  I,  that  the  down wash  is  finite  from 

a finite  doublet  so  that  a single  distributed  doublet  will  give  zero.  This  is 

not  the  case  for  but  while  this  derivative  has  infinities  at  the  Mach  Cone, 

dy 

the  integrated  effect  of  these  is  zero. 

In  order  to  obtain  the  pressure  distribution  LOO  , we  must  fill  the 
surface  of  the  wing  with  lines  of  doublets  at  various  angles  tu"  At  from 

d\a--c  to  h1 where  G is  the  tangent  of  the  half  angle  of  the  wdng.  The 

integral  equation  is  * 


III. 14 


, 

by  _ _ V X*~( M -Q  [ LO>)  dK 

' airpu  J.i  (h  x -S  )* 


¥e  interpret  III. 14  in  the  same  way  we  consider  a vortex  line  in  subsonic 
flow  because  of  the  similarities  of  the  infinite  velocities  occurin  in  the  two 
types  of  flows.  For  example  a line  vortex  in  subsonic  flow  located  at  )C  gives 

a c'ownwash. 


III. 15 


W- 


p 

Jltt  Ox,-*) 


<wd  *8/  d — 

dU  *TT(x,-*r 
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The  in  III. 14  due  a single  line  of  doublets  must  therefore  exhibit  an 
infinite  velocity  of  the  form 

To  solve  the  integral  equation  III. 14  i.e.,  we  Reed  the  function  LJ^which 
will  give  for  ||<Jc{and  which  is  symmetrical  with  respect  to  the  center- 

line,  we  tm«  L(K)=:  -A- . For  other  values  of  , we  can  utilize  the  results 

of  the  lift,ing-line  integral  eauation  of  Reference  3. 

In  order  to  evaluate  K j we  calculate  for  m>c  and  integrate 

O w- 

from  the  Mach  line  to  the  wing  through  the  singularity  at  the  wing  tip.  Performing 
this  integration  yields: 


III* 16 


aw  _ k ^ (y/x) 

iX  ~ JpU  X x C*)'*7* 


This  integration,  with  respect  to  X for^,  can  be  most  easily  performed 


by  letting  ^ to  yield: 

X 


. k C'  \ 

m.„  W-»«  = -5js£— JXX^pir 

/km 

We  integrate  III. 17  by  parts  to  obtain: 


III. 18 


voi  - ~ k L 1 ~ ✓ k«**m)Ava  . > rVaA i 


The  first  integral  is  zero  on  the  Mach  cone  and  infinite  at  X-C  • 

This  infinity  however  can  be  set  equal  to  zero  since  w must  be  Uo6  on  the 
airfoil  according  to  the  boundary  condition.  ^ explanation  of  uhis  singularity  is 
given  in  the  discussion  of  leading  edge  suction  contained  in  this  report. 
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Let  us  now  consider  the  second  inte  ral  and  let  r*  to  yielu: 


III. 19 


V lot,  - 


_ K f'k 


j±t_ 


This  integral  is  '-novTi  from  previous  work  to  be  an  elliptic  in>  e 
the  second  kind,  Reference  U,  P.  135.  Thus  -e  obtain: 


III. 20 


K ~ 


acxuV  p 


and  the  lift  distribution  is  even  by  the  ennression: 


III. 21 


LU0  = 


- ac iwVp. 


/c1-  K1  E(h-  cm’-Oc*-) 


This  well-known  result  is  contained  ir.  many  previous  references. 


ral  of 
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NON— STEADY  FhloSUKE  JGUELETS 


The  non-steady  pressure  doublet  can  be , derived  from  the  non-steady  pressure 
source  by  differentiation.  The  expression  for  the  non-steady  pressure  source 
is  given  in  reference  1 as  follows: 


III. 22 


|p  (X) 


i £ f(T*) 

47T  /x^  (Mx-»Xyx+Z>) 


where  € ftfj  represents  the  strength  of  the  source  as  a function  of  ^T^which  is 
related  to  and  ^ by  the  following : 

III.23  T = t MlC  -t-  V X-  (***  zjj 


C (M^i) 


Tliis  expression  means  that  each  point  of  the  fluid  X,  4/2.  receives  pressure 
impulses  at  time  ~fc  generated  by  the  sources  at  time T* . In  the  supersonic  case 
there  are  two  pulses  received  simultaneously,  a backward  moving  pulse  corresponding 
to  the  larger  value  of  and  a forward  moving  pulse  corresponding  to  the  smaller 
value  of 


For  an  oscillating  source  represented  by 
expression  yields: 


6 e 


the  pressure  source 


III.24 


(Mv,XvjVz?)i  C 


The  pressure  along  the  Mach  Wave  from  both  pulses  has  the  phase  lag  of 
so  that  wave  front  has  a sinusoidal  spatial  distribution  of  pressure  as 
might  be  expected. 
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The  pressure  doublet  is  iven  by  differentiation  vith  respect  to  21  as 


mentioned  above,  and  involves  four  terms  i.e., 


III.25 


- y d liSX  JLA1  V C-kJ  lUUi  1/^iUO  • to  • J 

{ 


li  -4-iii- 


Ul  } 

cw-D 


[x ^ c 06^+z1)) c ( w 

where  the  sign  means  the  sum  of  the  plus  ana  minus  terms. 

The  vertical  acceleration  of  the  air  particle  caused  by  the  doublet  is  given 
by  Newton’s  Law  . The  expression  for  |P  at  Z?<*is  the  same  as 

equation  III. 25  except  that  we  must  divide  by  and  set  21-0. 

Suppose  ve  desire  to  compute  the  W velocity  of  a particle  at  Xj9,2L“°  at 
some  value  of  time  "fc  . We  can  first  consider  the  increment  of  downwash  when  the 
particle  is  at  X^^ifOof  the  moving  system.  The  vertical  acceleration  is  given 
by  the  pressure  gradient  at  X|  but  with  a time  ~t ~~ increment  of  velocity 
is  os  that  the  expression  for  W is: 

-MV  ^ 

.26  W(x, 1-0 ,t)  = - -f|  ( x,  y,  2-*°,  t ^’)  & 


III 


were  X nei'ere  to  the  value  of  X at  the  Mach  cone.  The  final  formula  for  the 


downwash  caused  by  a doublet  is  then: 


III. 27 


x* 


18 
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The  solution  for  the  above  integral  has  not  been  obtained,  but  some  idea 
of  the  functions  involved  can  be  obtained  from  the  following  ap:  roach  which 
involves  an  integration  with  respect  toV|  for  the  source  before  developing  the 
doublet . 

In  this  approach  we  will  integrate  equation  III .24  with  respect  to 
between  the  two  values  of  ^ which  make  the  denominator  vanish.  In  reference  1, 
we  use  a substitution  similar  to  the  following: 


(MM) 


so  that 


III. 28 


coS  6 t ’ r*- 

\J  Sine 

dvi  * 

J <1 M1-! 


, cose 

tjrfjs  WL«  "•*■> 


-h  e 


Air  (/-V 

de 


] 


which  yields 


_ £ G 


m-29  1*f-'***>a 


i -<sy 


V CCMV,)  / 


Tire  pressure  caused  by  a line  of  doublets  is  obtained  from  III. 29  by 
differentiation  ie: 


III. 30 
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€ e 


2-  / W - I 


pp-  yj, 
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The  pressure  gradient  atZ-Ocan  be  obtained  iron:  III. 30  by  dividing  by  Z. 
and  setting  2.'°  • Using  Newton's  Law  as  explained  previously  for  the  oscill- 
ating down1,  ash,  we  obtain  the  following  eruption: 


j«(i  - %) 
W =--e  1 e 

3(30  1 MV  i 


* 


The  function  in  the  above  equation  could  be  tabulated  to  obtain  the 
oscillating  downwash  behind  a two— dimensional  lifting  line.  For  finite  lifting- 
lines a similar  expression  is  obtained  except  that  J0is  an  incomplete  integral. 
This  relation  shows  that  type  of  integrals  wuich  may  e expected  from  eouation 

III. 27* 
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LEADING  EDGE  SUCTION 

A suction  force  is  provided  when  air  flows  around  a corner  without  separation 
Suction  occurs  at  the  leading  edge  of  a subsonic  airfoil  or  the  leading  edge  of  a 
supersonic  airfoil  when  the  edge  is  swept  behind  the  Mach  cone  and  the  airfoil 
develops  lift.  A physical  picture  of  the  suction  force  is  important  in  the  design 
of  airfoils  for  subsonic  and  supersonic  flight. 

In  this  part  of  the  report,  a simple  mathematical  calculation  of  the  suction 
force  shall  be  given  based  on  linear  theory.  The  pressure  at  X^Zpaused  by  a 
subsonic  doublet  at  is  as  follows: 


III.32 


Q-fVQZ. 
4tr  *+  o-M*Xh-»r 


To  obtain  the  pressure  due  to  a line  of  doublets  extending  from  h~-i/o  to  n 
we  let  X-CS}'*)  ^ M -I  and  obtain 


III. 33 


which  gives: 


III. 34 


K,  — . L [ _dX 

« 4TT  J + 


3/a. 


whence : 


and 


> = 

"Ka 

- 


-u(y 

JTr 

_ L f x^- 

. L I^P-  3 ZX 

2TT 
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Using  Newton’s  Law 
the  following: 


2% 


and 


r-  we  obtain  for  U and  W , 


III. 35 


_ L<1$  j f-M'a- 


airpu  (x*4-(  l- 


111.36 


u - 


Ldi  /"hS ** 


140-M^IZO 


The  .elation  is  seen  to  hold  between  III. 34  and  III. 36. 

The  above  relations  for  Wand  U are  the  same  as  these  for  a two-dimensions  1 
vortex  in  compressible  flow  where in  III. 35  refers  to  the  downward  velocity  so 
that  sign  was  changed  in  the  process  of  going  from  III. 34  to  III  .35. 

For  a flat  plate  at  angle  o(  , we  use  III. 35  with  ©and  ^substituted  for  >c 
yielding  the  well-known  integral  eou?tion: 


III. 37 


W-  0.  - C-  L<t* 

Jb  (PTT 


From  reference  3 


ls  given  as  J3 J3  M v.-he^e  X=-bt*SO. 

In  order  to  calculate  the  suction  orce  we  enclose  the  leading  edge  with  a 


c'.irve  given  by  the  eouations 


* 


•S  — . On  t .is  curve  we 

b Vh*** 

calculate  Kf  and  W and  use  the  momentum  eouations  to  obtain  the  force  caused 


by  the  leading  edge.  In  this  tre  tment  §is  small  so  ti.au  the  range  of  integration 
of  0can  be  small  as  long  as  0 , Ti.e  expression  for  w and  V on  this  curve 

are  as  follows: 
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f3  -IS 

JL  r 


-[S  Cose  4 de 


III.38  W — A't*  

Tr  J0  [O1/!  4 ScoS^]1-*-  S*  W* 


III. 39 


so  that  III .40 


y . au^  f e s /y  i-fade 

^ Jo  +ScoS(f  J1  4 S^‘&lna'<ji> 

w-xffs-.„.  -W  .4;^» 


III. 40 


W3  _ au.et  CoS  4/x. 

KTT 


W / ina  mbMbMMi 

ra|  yr^ 


In  order  to  calculate  the  suction  force  from  III. 40,  we  need  the  pressure 
and  momentum  eouations  for  linearized  flow  retaining  terms  up  to  the  second 
order.  These  equations  are  derived  as  follows: 


III. 41  — » 


— y dv  -H^o  *vi ) du 


— c i — Uoul  \ 

where  p according  to  the  linear  theory  is  ^ — (*o  \ * 'ijfo*  ) 
Substituting  the  expression  for  p an  integration  of  III. 41  yields: 


III. 4-?  • 


- ±.  = UcU  + <£(>-« o%)  + 

Po 


For  the  momentum  equation,  we  oi.  tain: 


III. 43 


4^  =■  JplUof''Hz.  -(-^poVc^ 
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wliici;  alter  subsliLutiun  for 


i g la  z J 


111,44  4 "Tt*  3 -V  ju*<Jz  + 

The  suction  force  with  zero  pressure  round  t.  control  is 

m.45  £ - - - £ nvol% 

- 3 TT  bO,N;*- 

y hi1 


The  momentum  caused  by  the  pressure  is: 


111.46 


. (■* ▼ r+ir 

- * \ bit^\  (O 
at 


Therefore  the  expression  for  she  dry  force  comes  out  would  be  expected 
from  the  Imown  fact  that  no  pressure  rra£  exists  at  s i u..ni::  cr  eeds, 

hre  will  now  consider  the  suction  force  o a supersonic  lei  La  ving,  Ti;e 
lift  distribution  is  given  by  formula  III .21  . . o t.xit  U is  obtain^.*  by 

dividing  this  formula  by  — . 2 f>0^ 


III. 47 


VJ=  + 

v'te Wf  E ( /i- (mV, icd 


upper  surface 
lover  surface 
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Vmay  be  obtained 
Newton's  Law  erua valent 
written  as: 


integration  of  the  relation 


p*SS  'I?  • 


^ a#  ’ 

The  expression for  dV 

dx. 


or  its 
may  be 


III. 43 


av  + £**4*  %“*-  2l  «fe  * 1* 


Now  V-°^  the  leading  edge  at  L20  since  U=<*0  at  the 
leading  edge.  Integration  of  III. 43  from  0^  to yiel  .s: 


i.1. 49 


v — m.  Uo°<  ^/k 


The  singular  part  of  tins  integral  is  neglected  because  V goes  to  infinity 
li^-e  VI  at  the  leading  edge.  ontioncd  above  . 

At  the  leading  edge  ecuals  t,  and  ^?Cso  that  t,.e  resultant  velocity 

is  normal  to  the  leading  edge  and  has  a constant  value.  This  means  tltat  the  flow 
is  two-dimensional  in  a.  direction  normal  to  the  leading  edge.  The  value  of  this 
resultant  velocity  near  the  leading  edge  is 


III.  50  fuMvV 

or  letting  M ~ C-  & 

% **  ® 


Uoo*  c.  V ^ 


l/uMY*  =* 


\/c  V cx+ 1 

Wt  c £ ) 


Comparing  with  III. 40,  V and  w Snould  be  iven  by  the  relations 
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III. 51 


y-  1 Sm 

\Ti|  E 


.s  Uo^  _ jTc~  J^EFL-  Cas 

Voi  E 

where  the  Mach  number  normal  to  the  leading  edge. 

A verification  of  this  formula  at  ^>-0  can  be  obtained  by  considering  the 
singular  term  in  equation  111.18  letting'  v/jc^c  + S 


m„52 


w <*  ij±±1?LIE^S: 

E \Ti? 


{ 1- 

-t  MMII  ■»—  ■— 

/77c*~ 

The  suction  force  can  be  obtained  in  the  same  way  as  in  equation  111.45 


and  noting  that  ( ( — “ 


rr  ITpSA  or  the  total  suction  force  considering  the 

forward  component  is  for  both  sides: 


III. 53  f 


IT  9 
2. 


f Uo1^1  cx  yT*<> 

v <EV 


where  ]£  is  the  length  of  the  maximum  chord. 
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CONCLUSIONS 


1)  This  report  contains  the  derivation  of  the  doublet  flow  on  the  basis  of 
physical  reasohingi 

2)  The  physical  picture  of  doublet  flows  permits  one  to  comprehend  the 
sigularities  involved  and  to  better  appreciate  some  of  the  approximations  of 
linear  theory, 

RECCK^NDATICNS 

1)  Numerical  work  should  be  carried  out  to  obtain  the  effects  cf  unsteady 
motions  on  supersonic  dcwnwash. 

2)  Further  work  should  be  done  to  clarify  the  mechanism  involved  in  passing 
from  the  fixed  to  the  moving  doublet. 
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